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The Josephson effect [lj represents perhaps the 
prototype of macroscopic phase coherence and is 
at the basis of the most widespread interferom- 
eter, i.e., the superconducting quantum interfer- 
ence device (SQUID) [2j. Yet, in analogy to elec- 
tric interference, Maki and Griffin |3j predicted 
in 1965 that thermal current flowing through a 
temperature-biased Josephson tunnel junction is 
a stationary periodic function of the quantum 
phase difference between the superconductors. 
The interplay between quasiparticles and Cooper 
pairs condensate is at the origin of such phase- 
dependent heat current, and is unique to Joseph- 
son junctions. In this scenario, a temperature- 
biased SQUID would allow heat currents to in- 
terfere [4, 5j thus implementing the thermal ver- 
sion of the electric Josephson interferometer. The 
dissipative character of heat flux makes this co- 
herent phenomenon not less extraordinary than 
its electric (non-dissipative) counterpart. Albeit 
weird, this striking effect has never been demon- 
strated so far. Here we report the first experi- 
mental realization of a heat interferometer. We 
investigate heat exchange between two normal 
metal electrodes kept at different temperatures 
and tunnel-coupled to each other through a ther- 
mal 'modulator' [5j in the form of a DC-SQUID. 
Heat transport in the system is found to be phase 
dependent, in agreement with the original predic- 
tion. With our design the Josephson heat inter- 
ferometer yields magnetic-flux-dependent tem- 
perature oscillations of amplitude up to ~ 21 
mK, and provides a flux-to-temperature trans- 
fer coefficient exceeding ~ 60mK/$o a t 235 mK 
($ - 2x 10 15 Wb is the flux quantum). Besides 
offering remarkable insight into thermal trans- 
port in Josephson junctions, our results represent 
a significant step toward phase-coherent master- 
ing of heat in solid-state nanocircuits, and pave 
the way to the design of novel-concept coherent 
caloritronic devices. 

To figure out a Josephson heat interferometer we con- 
sider a symmetric DC-SQUID (i.e., a superconducting 
loop comprising two equal Josephson tunnel junctions 
with resistance Rj) composed of two identical supercon- 
ductors Si and S2 in thermal equilibrium and residing 
at temperatures T\ and T 2 , respectively (see Fig. la). 
For definiteness, we assume T\ > T 2 so that the SQUID 
is temperature-biased only. Within this assumptions, the 
total heat flow Q squid from Si to S2 becomes stationary 



and is given by [3J |3HZ] 

Qsquid{&) = 2Q qP — 2Qi 



cos 



7T$ 



(1) 



where the factor 2 accounts for two identical SQUID 
junctions, and <1> is the applied magnetic flux threading 
the loop. ^-dependence appears in Eq. only 
through the cosine term so that Q squid consists of a 
$o _ periodic function superimposed on top of a magnetic 
flux-independent component. In the above expression, 
Q qp (Ti , T 2 ) = ^ / °° desAft (e, T t )M 2 (e, T 2 ) [h (e, Ti) - 
/2(e, T2)] is the usual quasiparticles heat 
current [3J [5], whereas Qint{Ti,T 2 ) = 
5 3 | 7 / oo d ee ^ 1 ( e ,Ti)X2(£,r 2 )[/i( e ,Ti) h{e,T 2 )\ 
is the power flow due to interference between quasi- 
particles and the Cooper pairs condensate [3j IH [6j- 
|S]. fi(s,Ti) = [1 + exp^/fcsTj)]- 1 is the Fermi- 
Dirac distribution at temperature Tj (i = 1,2), 
M l (e,T l ) = H/V^-A^)^ 2 - A,(T0 2 ] is the 
BCS normalized density of states in the superconductors 
[ID], Mi(e,Ti) = A l (r 4 )/^ 2 -A. i (T l ) 2 6[e 2 - A^T,) 2 ] 
[3j HI [6HS] , Aj(Tj) is the temperature-dependent energy 
gap [TU], G(x) is the Heaviside function, fee is the Boltz- 
mann constant, and e is the electron charge. We note 
that both Q qp and Qmt vanish for T\ — T 2 , whereas Qi n t 
also vanishes when at least one of the superconductors 
is in the normal state. 

The implementation of our heat interferometer is 
shown in Fig. lb. The structure has been fabricated 
with electron-beam lithography and three-angle shadow- 
mask evaporation. It consists of a source and drain cop- 
per (Cu) electrodes tunnel-coupled to a superconduct- 
ing aluminum (Al) island (Si) defining one branch of 
a DC-SQUID. Source and drain junctions normal-state 
resistances are R s ~ 1.5 kO and Rd — lkfi, respec- 
tively, whereas the resistance of each SQUID junction is 
Rj ~ 1.3 kfi. Si is also contacted by an extra Al probe 
(S3) via a tunnel junction with normal-state resistance 
R p ~ 0.55 k£7, enabling independent characterization of 
the SQUID. Both source and drain are tunnel-coupled to 
a few external Al probes (vertical wires in Fig. lb) so 
to realize normal metal-insulator-superconductor (NIS) 
junctions, with normal-state resistance ~ 25 each, 
which allow Joule heating and thermometry [}J]. 

Below the critical temperature of Al (~ 1.4 K) Joseph- 
son coupling allows dissipationless charge transport 
through the SQUID. The SQUID voltage-current charac- 
teristics at 240 mK for two representative magnetic-flux 
values arc shown in Fig. lc. In particular, a well-defined 
Josephson current with maximum amplitude of ~ 226 



2 



,--? sou 



,(*) 



q2 



Eo 



oyo 








1/2 


ft 


t 


t 


ifl 




1 i 1 i 1 



thermometer 




FIG. 1. Josephson heat interferometer, a, Scheme of 
the device core: two identical superconductors Si and S2 kept 
at different temperatures Ti and T2 (with Ti > T2), respec- 
tively, are tunnel coupled so to implement a DC-SQUID. Rj is 
the normal-state resistance of each junction, $ is the applied 
magnetic flux threading the loop, whereas Qsquid{&) is the 
heat current flowing from hotter to colder superconductor, b, 
Scanning electron micrograph of the heat interferometer. Cu 
source and drain electrodes share the contact through K\O x 
tunnel barriers with an Al island (Si) defining one branch of 
a DC-SQUID. The other branch of the SQUID (S3) extends 
into a large volume lead to insure proper thermalization of its 
quasiparticles at the bath temperature Tt, a th- SQUID junc- 
tions are marked by crosses. Si is also in contact with an 
Al tunnel probe (S3) enabling independent characterization 
of the SQUID. NIS junctions made of Al are used as heaters 
and thermometers. See text for details, c, SQUID voltage 
(V) versus current (/) characteristics at two representative 
values of the applied flux. $0 corresponds to an applied field 
B = $0 1 A ~ 1 Oe, where A ~ 19.6 (an 2 is the ring area, d "In- 
dependent pattern of the SQUID critical current I c . Dashed 
line is the theoretical result for a DC-SQUID assuming ~ 0.3% 
asymmetry between the critical currents of the two junctions. 
Data in c and d were taken at 240 mK. d, Source and drain 
NIS thermometers calibration curves. Symbols indicate the 
measured voltage drop Vth at a bias current of 70 pA (source) 
and 40 pA (drain) versus Ti, a th- Dashed lines are theoretical 
results for a NIS junction. 



nA is observed at $ = 0. The magnetic-flux pattern of 
the SQUID critical current I c along with the theoreti- 



cal prediction [TU] is displayed in Fig. Id, and shows a 
nearly-complete supercurrent modulation which confirms 
the good symmetry of the SQUID. 

Thermal transport and, therefore, heat interference in 
the structure, arises from heating electrons in the source 
above lattice temperature (Tbath) so to elevate the quasi- 
particles temperature in Si (Ti) and create a temperature 
gradient across the SQUID. This hypothesis is expected 
to hold as the second branch of the SQUID (S2) extends 
into a large volume lead providing efficient thermalization 
of its quasiparticles at T oa th- Q squid will thus manifest 
itself leading to a $o-periodic modulation of drain elec- 
tron temperature {Tdrain)- 

Investigation of heat transport in our system is per- 
formed as follows. One pair of NIS junctions in the source 
is operated as a heater whereas a second pair is used 
to measure electron temperature {T source ) by applying a 
small DC bias current and recording the corresponding 
temperature-dependent voltage drop Vth (H E] ■ Analo- 
gously, another pair of NIS junctions is used to perform 
thermometry in the drain (see Fig. lb). Thermometers 
bias currents were optimized to achieve high sensitivity 
while limiting the impact of self-heating or self-cooling 
[S]. Figure le displays the calibration curves of source 
and drain thermometers versus T oat h, obtained by slowly 
sweeping the cryostat temperature from 235 mK to 750 
mK. The corresponding theoretical results for a NIS junc- 
tion [TU] are shown as well. 

Measurement of heat interference is done by stabilizing 
the cryostat temperature at a desired T\, a ih and heating 
the source up to a given T source . Tdrain is then recorded 
against a slowly-sweeping external magnetic flux. Figure 
2a shows Tdrain against <f> measured at 235 mK for in- 
creasing values of T source - Notably, Tdrain is $o-P er i°dic 
in $, as the Josephson critical current (see Fig. Id). 
As we shall argue, such a temperature modulation is of 
coherent nature, and stems from magnetic flux-control 
of Qsquid which is a hallmark of the Josephson effect. 
By raising T source leads to a monotonic enhancement of 
the average drain temperature over one flux quantum, 
(Tdrain), which follows from increased heat flow across 
the structure. On the other hand, the modulation ampli- 
tude STdraim defined as the difference between the maxi- 
mum and minimum values olTdraim turns out to initially 
increase and then tends to saturate at larger T source . In 
particular, STd ra in up to ~ 21 mK is observed corre- 
sponding to ~ 9% of relative modulation amplitude at 
235 mK. Tdrain modulation is observed also without in- 
tentional source heating, and might be related to a para- 
sitic power (~ 1— 5 fW) in the structure. The full T source - 
dcpcndcnce of (Tdrain) and STdrain are displayed in Fig. 
2b and confirm the above described behavior. We stress 
that heat interference manifests itself only owing to the 
existence of a finite temperature bias across the SQUID. 
Any voltage drop occurring at the Josephson junctions 
makes the phase-coherent component of Qsquid time- 



3 



dependent thus not contributing to time-averaged heat 
transport [3lll|8]. 

A relevant figure of merit of the heat interferometer 
is represented by the flux-to-temperature transfer coeffi- 
cient [5], T = dTdrain/d^, shown in Fig. 2c versus $ for 
a few selected T source . It turns out that |T| exceeding 
60mK/<i>o is obtained at 675 mK. Larger values might 
be obtained by lowering T oat h and by further optimizing 
the structure design. 

To account for our observations we have elaborated a 
thermal model sketched in Fig. 2d. We assume Si to ex- 
change heat at power Q source and Qdrain due to quasipar- 
ticle heat conduction with source and drain, respectively, 
at power Q squid with S2 and Q pro be with S3. Both S2 
and S3 are assumed to be thermalized at T&at/i- Further- 
more, drain electrons exchange energy at power Qdrain 
with Si, and at power Q ep with lattice phonons residing 
at Tbath [12] . The thermal steady-state of the system 
may be described by the energy-balance equations 
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; source 



QsQUId(&) + Qdrain = 

Qdrain ~t~ Qep O5 (2) 



where first equation accounts for thermal budget in Si, 
while the second one describes heat exchange in the 
drain. Ti and Tdrain can be determined under given con- 
ditions by numerically solving Eqs. ^ (see Methods 
Summary for further details). The model neglects heat 
exchange with photons due to mismatched impedance 
[T3HI7], electron-phonon coupling in Si owing to its re- 
duced volume and low experimental T oat h [18] . as well as 
phonon heat current [3]. 

Figure 2e shows the comparison between experiment 
and model displaying a few Td ra in{&) traces (solid lines) 
along with the theoretical behavior (dashed lines) . Anal- 
ogously, dashed lines in Fig. 2b show the predicted 
(Tdratn} and STdrain against Tsource- Although idealized 
the model provides reasonable agreement with our obser- 
vations, and grasps the relevant physical picture at the 
origin of heat interference in our system. 

The dependence on bath temperature is shown in Fig. 
3a (left panel) which displays Td ra in(&) at a few increas- 
ing T^h for Tgource set around 700 mK. Right panel of 
Fig. 3a shows the curves obtained from the model at 
the corresponding T oata . Besides leading to a monotonic 
enhancement of (Tdrain), by increasing T oat h yields sup- 
pression of STdrain and smearing of Td r ain(&) [5 which 
can be mainly ascribed to the enhancement of electron- 
phonon coupling in the drain [9] as well as to the influence 
of thermal broadening. STdrain ^2.5 mK is still observ- 
able at 450 mK, whereas the modulation disappears for 
Tbath ^ 500 mK. We emphasize that the latter is sub- 
stantially smaller than the temperature setting the disap- 
pearance of the Josephson effect in the SQUID (~ 1.4 K). 
Figures 3b-d show (Tdrain), STdrain, and the maximum 
of |T| versus T SOU rce, respectively, recorded at the same 
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FIG. 2. Behavior of the heat interferometer at 

235 mK. a, Flux modulation of Tdrain measured for several 
T source values, b, Modulation amplitude STdrain (left axis) 
and average temperature (Tdrain) (right axis) versus T source . 
Dashed lines are the results from thermal model (see discus- 
sion below and in the text), c, Flux-to-temperature trans- 
fer function T = dTd ra in/d<& versus <3> measured at a few 
selected Tsource- d, Idealized thermal diagram accounting 
for our setup. Si exchanges energy at power Q sour . ce and 
Qdrain due to quasiparticle heat conduction with source and 
drain, respectively, at power Q squid with S2 and Q pro be 
with S3. S2 and S3 are assumed to reside at T bat h- Drain 
electrons exchange energy at power Qdrain with Si, and at 
power Q e p with lattice phonons thermalized at bath tem- 
perature. Arrows indicate the direction of heat flows for 

Tbath < Tdrain < T\ < T SOU rce- e, Experimental T d rain{&) 

curves at a few selected values of T source along with the re- 
sults from thermal model (dashed lines). The full tempera- 
ture range in each panel is 28 mK, and the vertical division 
corresponds to 10 mK. 



Tbath as in panel 3a, along with the lines obtained from 
thermal model. Each of these quantities displays simi- 
lar qualitative behavior at different bath temperatures, 
and a smoothing of their characteristics is observed by 
increasing T oat h- This is consistent with the picture pro- 
vided by the model which captures the main features of 
the experimental data. 

Our results confirm what originally predicted almost 
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FIG. 3. Behavior of the heat interfer- 

ometer at different bath temperatures. a, 

Flux modulation of Tdrain recorded at different 
Tbath- From bottom to top, data correspond to 
T source = 675 mK, 700 mK, 690 mK, 700 mK, 700 mK. 
Left panel shows the experimental data, while the right one 
displays results from the thermal model. A sizable tem- 
perature modulation is still observable at 450 mK, whereas 
STdrain vanishes for Tbath > 500 mK. b, Average temperature 
(Tdrain) versus T SOU rce- c, Modulation amplitude STdrain 
versus T aour ce- d, Maximum value of |T| versus T source . 
Data in panels b-d were measured at the same T aa th as in 
panel a. Solid lines correspond to the thermal model. 



fifty years ago [3], i.e., the existence of a phase-dependent 
component in the heat current flowing through a 
temperature-biased Josephson tunnel junction. Besides 
offering insight into thermal transport in Josephson 
weak-links, this effect could represent a valuable tool 
toward phase-coherent manipulation of heat in solid- 
state nanocircuits [§1 ITM2"1"] . Yet, novel-concept coherent 
caloritronic devices such as heat transistors and ther- 
mal splitters could be envisioned which exploit phase- 
dependent heat transfer peculiar to the Josephson effect. 
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METHODS SUMMARY 

The structures have been fabricated with e-beam 
lithography and three-angle shadow-mask evaporation of 
metals onto an oxidized Si wafer through a suspended re- 
sist mask. In the e-beam evaporator, the chip is initially 
tilted at an angle of 28°, and 20 nm of Al are deposited 
to realize S2 and S3. The sample is then exposed to 380 
mTorr of O2 for 4.5 minutes to form the SQUID tunnel 
barriers after which it is tilted to 0° for the deposition 
of 25 nm of Al forming Si as well as heaters and ther- 
mometers probes. The chip is subsequently exposed to 
800 mTorr of O2 for 4.5 minutes to form heaters, ther- 
mometers, source and drain tunnel junctions. Finally, 30 
nm of Cu are deposited at 42° to realize source and drain. 

The magneto-electric characterization of the samples 
was performed down to 235 mK in a filtered 3 He cryostat. 
Current biasing of thermometers was obtained through 
battery-powered floating sources, whereas voltage and 
current were measured with room-temperature preampli- 
fiers. Flux-to-temperature transfer functions were mea- 
sured with low-frequency lock-in technique by superim- 
posing a small modulation to the applied magnetic field. 

In the energy-balance equations [see Eqs. 

§], Qprote = Q p q r P obe - QZf e EHS, Q p q r P obe = 

^ / °° deeM 1 {e)M s {e)[f l {e) - / 3 (e)], A/" 3 (e) = A/ 2 (e), 
M 3 (e) = M 2 {e), and / 3 (e) = f 2 {e). Further- 
more, Q source = ^jr J °° dzsM\ (s) [f s (e) - fl(s)], 
Qdratn = ^ J™ desMl (e) [/l (e) - f d (s)] , f s ( d )(e) = 
[l+exp(e /kBT source (drain))} , Qep = ^ (Tj ra m ~^bath) 

[12], S ~ 3 x 10 9 WK" 5 m- 3 is the electron-phonon 
coupling constant for Cu [9], and V ~ 2 x 10~ 20 m 3 is 
drain volume. For the numerical solution of Eqs. ([2]) 
we set the structure parameters as extracted from the 
experiment, and varied R p between ~ 100% and ~ 125% 
to match measured data. 
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